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ABSTRACT 



Context. With a Jupiter-mass planet orbiting at a distance of only 0.03 1 AU, the active K2 dwarf HD 189733 is a potential candidate 
in which to study the magnetospheric interactions of a cool star with its recently-discovered close-orbiting giant planet. 
Aims. We decided to explore the strength and topology of the large-scale magnetosphere of HD 189733, as a future benchmark for 
quantitative studies for models of the star/planet magnetic interactions. 

Methods. To this end, we used ESPaDOnS, the new generation spectropolarimeter at the Canada-France-Hawaii 3.6m telescope, to 
look for Zeeman circular polarisation signatures in the line profiles of HD 189733 in 2006 June and August. 

Results. Zeeman signatures in the line profiles of HD 189733 are clearly detected in all spectra, demonstrating that a field is indeed 
present at the surface of the star. The Zeeman signatures are not modulated with the planet's orbital period but apparently vary with 
the stellar rotation cycle. The reconstructed large-scale magnetic field, whose strength reaches a few tens of G, is significantly more 
complex than that of the Sun; it involves in particular a significant toroidal component and contributions from magnetic multipoles of 
order up to 5. The Ca ii H & K lines clearly feature core emission, whose intensity is apparently varying mostly with rotation phase. 
Our data suggest that the photosphere and magnetic field of HD 189733 are sheared by a significant amount of differential rotation. 
Conclusions. Our initial study confirms that HD 189733 is an optimal target for investigating activity enhancements induced by 
closely orbiting planets. More data are needed, densely covering both the orbital and rotation cycles, to investigate whether and how 
much the planet contributes to the overall activity level of HD 189733. 

Key words, stars: individual: HD189733 - planetary systems - techniques: spectropolarimetry - stars: magnetic field - stars: activity 



1. Introduction 

Planetary systems characterized by a giant planet at a few stel- 
lar radii from their parent stars (known as hot Jupiters) currently 
make up ~20% of all known extrasolar planets. Spectroscopic 
observations of stars hosting hot Jupiters recently revealed that 
chromospheric activity indices (such as core emission in Ca ii 
H&K and infrared triplet lines) are sometimes modulated with 
the orbital period of the giant planet (rather than with the rota- 
tion period of the star), suggesting that such giant p lanets may 
significantly boost the activ i ty leve l of the host star (iCuntz et al.l 
I2OOOI: IShkolnik et all l2003l l2005l). Ac cording to theoreticians 
(e.g., ICuntz et al.ll2000t lip et al.ll2004 . this modulation could 
result either from tidal effects (enhancing local dynamo action 
and activity at the surface of the star) or from magnetospheric 
interaction between the host star and its close-in planet (induc- 
ing reconnection events as the planet travels through the large 
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magnetic loops anchored in the stellar surface). If due to recon- 
nection, the resulting activity should depend mostly on the large- 
scale magnetic field of the host star, on the planet's magnetic 
field and on the orbital distance with respect to the Alfven radius 
of the host star; if due to tidal effects, the excess activity may in- 
dicate the presence of a "surface wave" of small-scale magnetic 
features comoving with the planet orbital motion and resulting 
from local turbulence and dynamo action enhancements induced 
by the planet pass-by. Photometric observations by the MOST 
space telescope of several hot- Jupiter systems suggest that stel- 
lar surface activity, in the form of cool spots comoving with 
the planet orbit al motion, co uld indeed be caused by the close- 
in giant planet d Walker et al. 2006); if confirmed, it would pro- 
vide additional evidence for magnetic interactions between giant 
planets and their host stars. Magnetic fields were also recently 
invoked as a possible formation mechanism for giant planets in 
close orbits around Sun-like stars (lKuruczll2007l) . 

Very recent observations of t Boo w ith the ESPaDOnS spec- 
tropolarimeter dPonati et aLlligQTllDon ati 2007), mounted at the 
Cassegrain focus of the 3.6m Canada-France-Hawaii Telescope 
(CFHT), revealed that surface magnetic fields are indeed present 
at the surface of planet-hosting stars and can be reliably detected 
through the Zeeman signatures they generate in spectral line pro- 
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files (ICatala et al.ll2007h . Despite a sparse temporal sampling of 
the rotational cycle, their data demonstrated the feasibility of 
measuring and mapping large-scale magnetic fields at the sur- 
faces of such stars, opening radically new and extremely promis- 
ing perspectives for observational studies of the magnetic inter- 
actions between active stars and their close-in giant planets. We 
therefore decided to initiate a new program aimed at detecting 
and modeling the magnetic topologies of stars hosting close-in 
giant planets, in order to provide theoretical models of planet- 
star magnetic interaction with quantitative constraints. With its 
short orbital period, the transiting planetary system HD 189733, 
discovered at Observatoire de Haute Provence through both pho- 
tometric and spectroscopic signatures (Bouchv et al,.2005.) . is an 
ideal candidate for our program. 

Most parameters of HD 189733 are now tightly constrained, 
thanks to successive very detailed spectroscopic and photomet- 
ric monitoring programs dBouchv et al.ll2005l:lBakos et al.H2006t 
IWinn et al.ll2007b a). The planet's mass and radius are equal to 
1 . 1 3 ± 0.03 Mjup and 1.156 + 0.046 Rjup respectively, while its or- 
bital period is 2.2185733+0.0000019 d. The system geometry is 
also well known, the inclination / of the orbital axis with respect 
to the line of sight being 85.76 + 0.29°. With a mass, effective 
temperature and radius respectively equal to 0.82 + 0.03 Mq, 
5050 + 50 K and 0.753 ± 0.023 Rq, the K2V parent star is 
slightly metal poor and active, with a line-of-sight projected 
equatorial rotation velocity vsin/ = 3.0 + 0.2 kms"'. In ad- 
dition to the photometric eclipses induced by the transiting ex- 
oplanet, HD 189733 exhibits intrinsic photometric variability; 
this is likely the result of its activity and the dark spots on 
its surface, modulating the amount of light received from the 
star as rotation carries spots in an d out o f the observe r 's view 
dHibrard & Lecavelier des Etangsl l2006t IWinn et alj l2007at 
lMatthews«2007i:iCroll e t al. 2007 ). From Hipparcos obseryations 
collected over 3 yr, [Hebrard & Lecaveher des Etangsl (l2006l) 
confirm that intrinsic variability is indeed the main source of 
out-of-eclipse photometric fluctuations on HD 189733, and that 
these fluctuations are modulated with an average rotation pe- 
riod of 11.8 d. From 93 nights of ground-based photometry 
secured from 2005 October to 2006 July, Winn et al. (2007a) 
show that the light curve of HD 189733 is evolving signifi- 
cantly on a timescale of 20 d. Using their last 30 d of obser- 
vations (showing the highest amplitude modulation), they derive 
a rotation pe riod of 13.4 ± 0.4 d, significantly longer than that 
of iHebrard & Lecavelier de s Etangs (2006). Nearly-continuous 
space-based photometry of HD 189733 was obtained over 21 d 
in July and Au gust 2006 w i th MOS T, only 50 d after the latest 
observations of IWinn et al.1 (l2007al) : from these data, the rota- 
tion period was found to be 11.8 d to within better than 0. 1 d 
dMatt hews 2007; Cr oll et al.ll2007l) . We suspect that the discrep- 
ancy between these different estimates is not due to measurement 
errors, but rather to temporal variations in the spot pattern cou- 
pled to a signi ficant amount of differential r otation over stellar 
latitudes (e.g., 'Hall''l99l|; ICatalaetZ 2007). We return to this 
issue later in the paper; in the meantime, we adopt an average 
rotation period of 1 1.8 d for HD 189733. Table[T]summarises all 
relevant system properties that are used in this paper and corre- 
sponding references. 

This paper presents the first spectropolarimetric study of 
HD 189733. Our aim is to detect and monitor the Zeeman sig- 
natures in the spectrum of the host star to map its large-scale 
magnetic topology and study the magnetic interaction between 
the star and its close-by giant planet. In Section |2] we describe 
the observations, while in Sect. [3] we detail the modeling of the 
spectropolarimetric data. We discuss the activity level of the star 



Table 1. Star and planet parameters (and associated error bars) 
for the transiting system HD 1 89733. Corresponding references 
are: " Bouchy et al. (2()05), iBakos et al.,.2006. ,Winn et all 
I2007bl "' ICroU et al.ll2007i 



Planet mass (Mjup) ^ 


1.13 (0.03) 


Planet radius (Riup) 


1.154 (0.032) 


Orbital period (d) ^ 


2.2185733 (0.0000019) 


Semi-major axis (AU) 


0.031 (0.001) 


Stellar mass (Mq) " 


0.82 (0.03) 


Stellar radius (Rq) ^ 


0.753 (0.02) 


Stellar effective temperature (K) " 


5050 (50) 


Stellar v sin / ( km s" ' ) ^ 


2.97 (0.22) 


Stellar rotation period (d) 


11.73 (0.07) 


Orbit inclination (°) * 


85.76 (0.29) 



and its possible variation with orbital and rotational phase in 
Section |4] In Sect. |5] we examine spectra taken during transit 
for potential spectroscopic signatures of the planet's atmosphere, 
and in Sect. |6] we derive upper limits on the detection of starlight 
reflected by the planet. Finally, in Sect. |7]we summarise the re- 
sults , discuss their implications and suggest directions to extend 
our study in the future. 

2. Observations 

Spectropolarimetric data of HD 189733 were collected with 
ESPaDOnS and CFHT, during 2 separate runs in 2006 June (10- 
13) and August (05-12). With ESPaDOnS, the complete optical 
spectrum (from 370 to 1000 nm) is collected in a single expo- 
sure at a spectral resolution of about 65,000; from sequence of 4 
subexposures taken in different configurations of the polarime- 
ter, one can obtain both the intensity and polarisation spectrum, 
in either circular (Stokes V) or linear (Stokes Q or U) polari- 
sation state (Donati et al. 1997; Donati 2007). Circular polari- 
sation spectra give access to the Stokes V Zeeman signatures 
and thus to the longitudinal field component averaged over the 
visible stellar hemisphere; they allow us to retrieve the parent 
large-scale magnetic field at the surface of the star (see Sec.[3]l. 
Linear (Stokes Q and U) polarisation can potentially inform us 
of the amount of stellar light reflected off the dusty planetary at- 
mosphere; this signal (whose typical size is of order a few tens 
of ppm with respect to the stellar flux and depends on the planet 
albedo) is expected to be modulated by the planet orbital phase 
(see Sec. 16]) and to reach its maximum size near the planet orbital 
quadrature. Eight Stokes V sequences at 3 main epochs, and 7 
Stokes V sequences at 6 main epochs were collected during the 
first and second run respectively, the 4 remaining sequences cor- 
responding to Stokes Q and U sequences. Signal to noise ratios 
(S/N) range from 600 to 1 100 per 2.6 km s"' velocity bin, with 
typical values between 900 and 1000. The complete journal of 
observations is given in Table |2l 

Orbital and rotational phases (listed in Table [2]l were com- 
puted respectivel y using the ephemerides (see lWinn et alj2007al ; 
ICroll et alJl2007l) : 

Tq = HJD 2, 453, 988.80336 + 2.2185733£orb (1) 

Tq = HJD 2, 453, 988.80336 + 1 1 .8£rot (2) 

A full planet transit was observed on June 13, with se- 
quences 10 and 11 corresponding to transit (ie orbital phase 
equal to 0.0+0.017), while sequences 9 and 12 correspond to 
phase ranges immediately preceding and following transit. 
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The raw data were reduced using Libre-ESpRIT, the ded- 
icated and fully automatic spectropolarimetric reduction pack- 
age available to ESPaDOnS users at CFHT (Donati et al. 1997; 
iDonatil 120071) . From collected calibration exposures and stel- 
lar frames, Libre-ESpRIT automatically extracts wavelength- 
calibrated intensity and polarisation spectra with associated er- 
ror bars at each wavelength pixel, shifted to the heliocentric 
rest frame and normalised to a unit continuum. In particular, all 
spectra are automatically corrected for spectral shifts resulting 
from instrumental effects (e.g. mechanical flexure, temperature 
or pressure variations) using the numerous telluric lines present 
in each spectra as a stable wavelength reference. Through cross- 
correlation, we estimate the relative shift of telluric lines with 
respect to the observer's frame; we assume that this shift is en- 
tirely attributable to instrumental effects and simply remove it 
by translating the spectra and recenter the telluric lines in the 
observer's rest frame. This is of course only an approximation; 
telluric lines are actually not a perfectly stable reference in the 
observer's rest frame, eg due to the effect of strong winds in 
the upper atmosphere (where telluric lines form) that can po- 
tentially shift these lines by as much as a few tens of ms"'. In 
practice, we find that our method yields a radial velocity dif- 
ferential rms precision of about 15 ms ' within a sin gle run, 
with s ystematic shifts of up to 30 m s ' from run to run jPonatil 
l2007h . This method is less accurate than the classical technique 
which consists in recording the spectrum of a reference ThAr 
lamp simultaneously (and interleaved) with that of the main star, 
yielding typical differential radial velocity accuracies of a few 
ms"'; our method is nonetheless accurate enough to detect the 
planet-induced radial velocity variations of HD 189733 (see be- 
low) and to ensure a reliable rotational modulation analysis of 
the detected Zeeman signatures (see Sec.|3]l. 

Least-Squares Deconvolution (LSD, iDonati et al.ll99 7*) was 
applied to all spectra, using a line pattern derived from a Kurucz 
model atmosphere with solar abundances, and effective temper- 
ature and logarithmic gravity set to 5000 K and 4.5 respectively. 
This line pattern includes most moderate to strong lines present 
in the optical domain (those exhibiting central depths larger than 
40% of the local continuum, before any macroturbulent or ro- 
tational broadening, about 7,000 throughout the whole spectral 
range) but excludes the very strongest, broadest features such as 
the Balmer lines, whose Zeeman signature is strongly smeared 
out compared to those of narrow lines. The output of LSD for 
each sequence of polarisation spectra is a pair of average line 
profiles: either intensity (Stokes /) and circular (Stokes V) polar- 
isation profiles, or linear (Stokes Q and U) polarisation profiles, 
depending on the type of observation sequence. We will refer to 
these as LSD profiles throughout the paper. The typical multi- 
plex gain in S/N for polarisation profiles is 30, implying noise 
levels in LSD polarisation profiles as low as 30 to 40 ppm (ie 
3 to 4 X 10"^ in units of the unpolarised continuum) depending 
essentially on the input spectrum quality. An example of such 
profiles is shown in Fig.|2] 

Stokes V Zeeman signatures are detected for all circular po- 
larisation observations (15 profiles altogether). From the first 
moment of the Stokes V LSD profiles, we derived estimates 
of the line-of-sight component of the magnetic vector averaged 
over the visi ble stellar hemisphere, ie the so-called longitudinal 
field Be (Donati et alj[l997) . The values obtained are listed in 
Table|2]and typically range from -8 to -1-8 G. Note that longitu- 
dinal field values only represent a small fraction of the informa- 
tion stored in Stokes V profiles. For example, while the Stokes V 
signature corresponding to sequence 16 is clearly detected, fea- 
turing a signal with a full amplitude of 0.06% (ie 20 times larger 



HD 189733, LSD profiles, 2006 Jun 11 
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Fig. 1. LSD unpolarised (bottom) and circularly polarised (top) 
profiles of HD 189733 on 2006 June 11. AU LSD profiles ob- 
tained on this night (corresponding to observation sequences 3 
to 5, see Table|2]) were averaged together, further decreasing the 
relative noise level (in the polarisation profile) down to 19 ppm. 
In this particular example, the full amplitude of the Stokes V 
signal is 0.12%, ie 60 times greater than the rms noise level. A 
Zeeman signature is obviously detected, corresponding to a lon- 
gitudinal field value of about 8 G. Note that the Stokes V profile 
was expanded by 100 and shifted upwards for graphics purposes. 



than the corresponding error bar), the corresponding longitudi- 
nal field value is consistent with G within the noise level and is 
thus unable to diagnose properly the field detection we obtained. 
Nonetheless, longitudinal-field values are convenient summary 
statistics in several respects, and in particular can provide useful 
estimates of the average magnetic flux over the stellar surface 
when the field polarity is roughly uniform over the visible stellar 
hemisphere. In particular, these values indicate straightforwardly 
(and not unexpectedly) that Be correlates very poorly with orbital 
phase. Sequences 3 and 12, for example, correspond to very sim- 
ilar orbital configurations, but yield widely discrepant Be values. 
On the other hand, Br varies rather smoothly with stellar rotation 
phase, changing from positive to negative values throughout the 
first observing run. A more complete modeling of the magnetic 
topology from our collection of LSD Stokes V profiles is pre- 
sented in Sect. [3] 

From Gaussian fits to each LSD Stokes / profile, we derive 
radial velocities of HD 1 89733 throughout our 2 runs. The values 
we obtain (listed in Table |2]i are in g ood agreement with expec- 
tations (using the orbital solution of Bouchv et al.ll2005l and the 
ephemeris of Winn et al. 2007a given in eqn.[T) once the data are 
globally shifted by -0.13 kms"' and -0.23 kms"' for the June 
and August sets respectively (see plot on Fig.|2]i. Residuals to the 
model are equal to 15 and 10 m s ' rms for the June and August 
run respectively, confirming that this is indeed the typical radial 
velocity differential accuracy that ESPaDOnS can reach. While 
the global shift of abou t 0.2 kms~' bet ween the ESPaDOnS and 
the ELODIE data (.Bouchv et al.1 12005 ) may be attributable to 
differences in the radial velocity calibration procedures between 
the 2 instruments, the 100 m s"' relative shift between our June 
and August data sets is significantly larger than the upper limit 
corresponding to the long-term stability of ESPaDOnS (about 
30 m s"'). It suggests that HD 189733 could host more than one 
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HD 189733, ESPaDOnS data, 2006 Jun & Aug 
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Fig, 2, Radial velocities of HD 189733 derived from ESPaDOnS 
spectra collected in 2006 June (crosses) and August (squares) as 
a function of orbital ph ase (using the ephern eris of eqn.[T]i. The 
radial velocity model of lBouchv et alj (l2005l) is also plotted (full 
line). Note that the measurements were shifted by -0.13 km s"' 
and -0.23 kms ' for the June and August data respectively, to 
compensate for a global shift with respect to the model predic- 
tions (see text). Individual 1 - cr error bars of 15 ms"' are also 
plotted. Note that the deviating point near phase 1 .0 is compat- 
ible with th e Rossiter-Mc Laughlin effect dBouchv et al.ll2005t 
IWinn et al.ll2007a) . 



giant planet and feature at least another one at a typical distance 
of 0.5 AU. This will be investigated through future long-term 
spectroscopic monitoring. 

Finally, we perform a bisector analysis from normalised LSD 
Stokes / profiles and measure the bisector span, from average 
velocities corresponding to 10 - 40 % and 55 - 85 % of the full 
LSD profile depth. The error bar on the bisector span is generally 
of the order of the error bar on the radial velocity estimate, about 
15 m/s in the present case. 



3. Modeling the magnetic topology of HD 189733 

To model the magnetic topology of HD 18 9 733, we use our 
magn etic-imaging code (Brown et al] 119911: Donati & Brownl 
Il997h in its most recent implementation dDonati et al.l l2006l) . 
While still based on maximum-entropy image reconstruction, 
this latest version reconstructs the field topology as a spherical- 
harmonic decomposition, rather than as a series of independent 
magnetic-image pixels. One obvious advantage of this method 
is that we can impose a priori constraints on the field topol- 
ogy - e.g., that the field is purely poloidal, or purely toroidal, 
or more generally, a combination of both. Another important ad- 
vantage of this formalism is that both simpl e and complex ma g- 
netic topologies can easily be reconstructed dDonati et alj|2006) . 
whereas the original method failed a t reconstructing sim ple mag- 
netic geometries (such as dipoles: [Brown et"anil99 1'). Details 
about how we effected this approach are given in|ponati et al. 
d2006l) . In this section, we assume that the Stokes V profiles are 
modulated with the rotation period of 11. 8 d; we further assume 



that the rotation spin axis of the star coincides with the orbital 
axis, inclined at 85.7° with respect to the line of sighfl 

To describe the LSD St okes V profil e s of H D 18 9733, we use 
the sam e simple model of iDonati et alj dl997h and lDonati et al.l 
d2006h . In this model, the local unpolarised line profile is de- 
scribed with a Gaussian, whose full-width-at-half-maximum is 
set to 7 kms"' and whose equivalent width matches that of the 
observed LSD Stokes / profiles (3.2 kms '); inte grating over 
the ste llar surface and taking vsin/ = 3.0 kms ' dWinn et al.l 
I2007al) enables us to reproduce the observed Stokes / profiles 
fairly well, once the model profiles are broadened with an in- 
strumental broadening of 4.6 kms"'. We then assume that the 
local Stokes V profile is proportional to the first derivative of 
the local Stokes / profile (weak field approximation), the scal- 
ing factor depending in particular on the local longitudinal field 
component and the Lande factor of the average lin e (set here to 
1.25). For more details, the reader should refer to iDonati et al.l 
dl997h . 

Given the fairly low v sin / of HD 189733, it is obvious that 
the spatial resolution we can achieve at the surface of the star 
is rather limited. From the width of the local line profile (about 
7 kms"'), we derive that there are about 3 resolution elements 
across the stellar equator, implying that using spherical harmon- 
ics expansions with ^ < 5 is sufficient for our needs. It also im- 
plies that most small-scale magnetic features potentially present 
at the surface of HD 189733, e.g., in the form of small bipolar 
spot groups like those peppering the surface of the Sun, will re- 
main totally undetected in spectropolarimetric data sets; this re- 
flects the fact that close-by opposite magnetic polarities mutually 
cancel their respective effect in the Zeeman polarised signatures, 
that are sensitive to the vector properies of the field. 

The maximum-entropy fits we obtain in this context for each 
of our 2 data sets are shown on Fig. [3] while the correspond- 
ing maps (assuming a magnetic field featuring both poloidal and 
toroidal components) are shown on Fig. [3] Both data sets are fit- 
ted down to the noise level and correspond to a reduced of 
order unity. Note that in both cases, the data set is only fragmen- 
tary, covering no more than 20% and 50% of the rotation cycle 
for the June and August sets respectively. It explains in partic- 
ular why only little magnetic flux is reconstructed at the stellar 
surface in regions not (or only marginally) constrained by obser- 
vations. This is due to the fact that the imaging code is tailored to 
produce the simplest surface magnetic topology compatible with 
the data; it therefore spontaneously biases the solution towards 
topologies featuring little magnetic flux on poorly observed re- 
gions of the stellar surface, provided that (i) this is compatible 
with observations at other rotational phases and that (ii) the as- 
sumed field topology is complex enough to allow a mixture of 
non-magnetic and magnetic regions (i.e., the spherical harmon- 
ics expansion describing the field topology includes orders sig- 
nificantly larger than 1). For this reason, we emphasise that the 
modeling we carry out in this paper is only preliminary and de- 
serves further confirmation and assessment through more com- 



' While lWinn et alj j2007al) find that t he spin and orbit axis are most 
likely aligned in HD 189733. ICroll et al.l ( l2007h conclude in favour of a 
moderate spin-orbit misalignment from their very accurate photometric 
data set. However, the latter mention that their result only holds if no 
differential rotation is present at the surface of the star, an assumption 
whose validity is fairly uncertain (see end of Sec. [3](. We thus assume 
here that the spin and orbit are aligned; a 20° change in the assumed 
inclination has however very limited impact on the derived magnetic 
topology. 
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Table 2. Journal of spectropolarimetric observations. The successive columns list the sequence number (col. 1), the UT date (col. 2), 
the UT Heliocentric Julian date (col. 3), the orbital phase (using ephemeris of eqn.[T] col. 4), the rotational phase (using ephemeris of 
eqn.|2] col. 5), the measured polarisation state (col. 6), the total exposure time (col. 7), the peak signal to noise ratio in the spectrum 
(per 2.6 kms"' velocity bin, col. 8), the longitudinal magnetic field B[ whenever Stokes V is measured (with corresponding Icr 
error bars, col. 9), the equivalent width S ca of the residual emission core in the Ca ii H&K line (col. 10, see Sect. |4]i and the radial 
velocity Vrad (col. 11). 



seq. 


date (UT) 


HJD 


0orb 


0rot 


Stokes 


exp. time 


peak S/N 


Be 


5ca 








(2 453 OOO+'l 








(sec) 


( npr niYpl^ 


CGI 


(kms"') 


(kms"') 


I 


7006 Tun 10 


897 01115 


0.625 


0.221 


n 


4x960 


980 




-0.56 


-2.105 


2 


2006 Tun 1 


897 05805 


0.647 


0.225 


u 


4x960 


1040 




-0.82 


-2.089 


3 


9006 Tnn 1 1 

£^\J\J\J J nil 1 1 


fiQl 95814 


0.052 


0.301 


Y 


4x900 


1080 


+7.2 + 0.6 


-0.1 1 


-2.292 


4 


9006 Tnn 1 1 

^KjyjKJ J nil 1 1 


898 00969 


0.072 


0.305 


Y 


4x900 


1000 


_|.7 7 + 7 


-0.53 


-2.321 


5 


9006 Tnn 1 1 

j—\J\J\J J nil 1 1 


898 04669 


0.092 


0.309 


Y 


4x900 


930 


+7 9 + § 


-0.16 


-2.338 


6 


2006 Jun 12 


898.98463 


0.515 


0.388 


V 


4x900 


830 


+2.4 ±0.8 


+0.20 


-2.191 


7 


2006 Jun 12 


899.03159 


0.536 


0.392 


Q 


4x1000 


1040 




+0.27 


-2.166 


8 


2006 Jun 12 


899.08036 


0.558 


0.396 


u 


4x900 


1140 




+0.11 


-2.154 


9 


2006 Jun 13 


899.99747 


0.971 
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Fig. 3. LSD Stokes V profiles of HD 189733 (thin line) and corresponding maximum-entropy fit (thick line) for our 2006 June (left 
panel) and August (right panel) data sets. The rotational phases of observations (as derived from the ephemeris of eqn.|2]) are noted 
to the right of each profile. A 3-cr error bar is also depicted to the left of each profile. 
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Fig. 4. Maximum-entropy reconstructions of the magnetic topology of HD 189733 as derived from our 2006 June (left panel) and 
August (right panel) data sets, assuming that the global field can be expressed as the sum of a poloidal field and a toroidal field. 
The three components of the field are displayed from top to bottom (flux values labelled in G). The star is shown in flattened polar 
projection down to latitudes of -30°, with the equator depicted as a bold circle and parallels as dashed circles. Radial ticks around 
each plot indicate phases of observations (as derived from the ephemeris of eqn. 2). 



plete data sets. However, a number of reliable conclusions can 
be derived already from the present data. 

The first point we can address concerns the topology of the 
field. By trying to fit the data with different field configurations 
(eg a poloidal field only, or a general poloidal plus toroidal field 
combination), one can check which topology is more likely to be 
present at the surface of the star In the present case, we find that 



the field of HD 189733 most probably features both a poloidal 
and a toroidal component at the surface of the star Although 
both options give equivalent likelihoods in terms of the quality 
of fit to the data, the inclusion of a toroidal field component pro- 
duces a magnetic configuration with significantly lower contrast 
(hence greater entropy and higher prior probability), yielding a 
higher posterior probability. This is particularly obvious for the 
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August data set, for which the poloidal plus toroidal field con- 
figuration that fits the data at the noise level is almost twice as 
weak on average as that obtained when assuming that the star 
hosts a purely poloidal field (not shown here). This conclusion 
can be infeiTed straightforwardly from the corresponding Stokes 
y data set (see Fig. [3]). At this epoch, the mean Stokes V profile 
averaged over the whole series, as well as more than half of the 
individual profiles (eg that at phase 0.620), are clearly more-or- 
less symmetric about the line center. This is the characteristic 
signa ture of an azimu thal field ring encircling most of the star 
(eg, iPetit et al.ll2005l) . The same conclusion may be deduced 
from the June data set, although the evidence is weaker than for 
the other (more complete) data set. In particular, the portions of 
the azimuthal field ring reconstructed at both epochs feature the 
same (ie clockwise or negative) polarity (see Fig. O, strength- 
ening the evidence that this toroidal field ring is indeed a real 
component of the magnetic topology of HD 189733. 

We can also be fairly confident that the magnetic field at the 
surface of the star is significantly more complex than that of the 
Sun. Magnetic field models with a spherical harmonic expan- 
sions restricted to ^ < 2 produce significantly worse fits to the 
August data set, with reduced values always larger than 2.5. 

The modulation of the Zeeman signature on the August run 
looks very different from the regular shift from one main po- 
larity to the other in about half a rotation period (as observed 
from most chemically peculiar stars. Wade et al. 2000a) that a 
simple dipolar-like topology would produce. Similarly, we ob- 
serve in June a complete transition from one polarity to another 
in no more than 15% of the rotation cycle, which again can- 
not be reconciled with a simple dipolar-like magnetic topology. 
This is also straightfowardly visible from the longitudinal field 
values that we obtain (see Table |2]i, whose temporal variation 
is in cons istent with expectations from a simple field topology 
(eg, IWade et al. 2000b). Although the data at both epochs are 
too sparse to give a definitive picture of the global magnetic 
topology, we can nevertheless infer that it includes several (pre- 
sumably radial) magnetic field spots of opposite polarities across 
the equator. The typical magnetic intensity within these spots is 
a few tens of G (and possibly up to 40 G in the azimuthal field 
ring), significantly greater than the large-scale field strength of 
the Sun (about 5 G at the pole). 

Further evidence of complexity comes from the apparent 
variability of the field topology on a time scale of a few weeks. 
This is directly visible on the reconstructed maps, where the 
magnetic spot distributions at both epochs share little in com- 
mon. Although this is not totally unexpected given the very in- 
complete phase coverage, this is also the case around phase 0.50, 
a portion of the stellar surface that was covered at both epochs. 
This is particularly obvious in the Zeeman signatures collected 
around phase 0.48 in 2006 June, which are obviously incompat- 
ible (both in shape and polarity) with that obtained at phase 0.47 
in 2006 August and with those obtained at nearby phases (0.38 
and 0.55). This is further confirmed by trying to fit simultane- 
ously both data sets with a single field topology. No acceptable 
fit to the data can be obtained for values of the period ranging 
from 11 to 14 d. Reduced values lower than 2 can only be 
obtained for some rotation periods between 12.2 and 13.4 d. For 
a period of 12.2 d, for instance, the Zeeman signatures at phases 
0.47 and 0.55 in the August data set are approximately matched 
with those at phases 0.30 and 0.38 in the July data set. We must 
therefore conclude that the magnetic topology evolved signifi- 
cantly between our two runs, e.g., under the effect of differential 
rotation. If we assume that the variability was not so large as to 
render the magnetic topologies at both epochs fully uncorrected. 



we find (by minimising the misfit between the two data sets) that 
the average recurrence period of the Zeeman signatures must be 
significantly longer than 1 1.8 d. 

This latter conclusion confirms reports in the literature from 
photometric monitoring that variations in light-curve morphol- 
ogy occur on HD 189733 on timescales of a few tens of days. 
As far as our data indicate, the recurrence period of the Zeeman 
signatures is longer than the 11.73-day period retrieved by 
Croll et al. (2007), but shorter than the 13.4-d period found by 
Winn et al., (,2007a ). We suspect that all these discrepant periods 
actually reflect the fact that HD 189733 is differentially rotating; 
the recurrence period derived from the various data sets traces 
essentially the rotation period at the latitudes of the main bright- 
ness/magnetic features that are present at the surface of the star at 
the time of the observations. For instance, all data sets could po- 
tentially be reconciled by assuming that the equator and pole of 
HD 189733 are rotating with periods of 1 1.8 and 13.5 d respec- 
tively, i.e., a latitudinal photospheric shear similar to that of the 
Sun. Providing this is indeed the correct interpretation, it would 
imply that differential rotation at the surface of HD 189733 is at 
least as strong as that of the Sun, i.e., with a lap time between 
the equator and pole of less than 100 d. Additional observations 
are obviously needed to confirm this option; spectropolarimetric 
data collected over a timescale of about 15 d could settle this 
issue fairly easily. 

4. Intrinsic stellar vs planet-induced activity 

HD 189733 is known as a v ariable star. I ts 5h k index value is 
0.525 from the catalog of Wrig ht et all 12004), implying that 
HD 189733 is among the 10% most active K dwarfs of the cat- 
alog, and among the 20% most active stars of the whole catalog 
(including 1,200 stars altogether). To investigate how the activity 
of HD 189733 varies with time, we searched the emission cores 
in the Ca ii H&K lines for potential variability. To maximise the 
accuracy, we computed a mean profile from both Ca ii H&K 
lines at each phase, and subtracted from each profile the aver- 
age over the whole data set. The resulting residuals are shown 
in Fig.|5] Variations of up to +10% of the core emission flux are 
observed throughout both runs. Changes in the amount of core 
emission are derived by measuring the equivalent widths of the 
residual profiles of Fig.|5] The coiTesponding values are listed in 
Table|2]and plotted in Fig.|6] When the H and K lines are treated 
separately, the same behaviour is observed in both cores, with 
similar amplitudes. 

The first obvious (and unsurprising) finding is that the activ- 
ity of HD 189733 (as witnessed through the Ca ii H&K emission 
cores) is intrinsically variable on time scales as short as 1 hr. On 
June 1 1 for instance, the activity noticeably increased between 
the first and second polarisation sequences, then dropped again 
on the third polarisation sequence. This sudden burst of activity 
is probably due to a small flare triggered by a local reconnection 
event; it is however impossible at this stage to speculate about 
the cause of this event given the very sparse data sampling that 
we have. 

The longer-term variations in core emission are apparently 
not related to the orbital phase, as one can also see fairly clearly 
from Table|2] While the activity is lower than average on June 1 1 
(just after the transit), it is much stronger than average about 2 
days later, both before, during and after the transit. More gen- 
erally, very little correlation is observed between the activity 
level and the orbital phase, for either data set. Activity corre- 
lates significantly better with stellar rotation phase (see Fig.|6]l, 
though intrinsic variability generates a significant dispersion of 
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data points. Figure [T] shows a fit to the nightly-averaged June 
data (top panel), describing the average rotational modulation in 
the activity level of HD 189733. Residuals with respect to this 
average rotational modulation are plotted against the planet or- 
bital phase in Fig. |7] (lower panels). These residuals both show 
a large dispersion of individual measurements and no obvious 
correlation with the planet orbital phase. The available phase 
sampling is however still unsufficient; we thus cannot conclude 
from our data whether the planet is responsible for activity en- 
hanc ement along the interaction models propos ed in the litera- 
ture jCranmer & Saaiil2 007; Mcl vor et al.ll2006l) . 

The variations in Ca ii H&K emission with stellar rota- 
tion phase appear to be in good mutual agreement for both 
runs(see Fig. |6l). Since the activity presumably relates to the 
large-scale magnetic topology, it suggests that the intrinsic vari- 
ability of the field structure as detected through the change in 
the Zeeman signatures (see Sect. [3]) is only moderate and local, 
leaving the global field mostly unaffected. It also suggests that 
the rotation period at the footpoints of the magnetic loops that 
confine and heat the chromospheric plasma to produce the ob- 
served Ca n H&K emission is close to 11.8 d, i.e . consistent 
with that derived from the MOST photometric data dWinn et al] 
l2007 a) but smaller than the recurrence timescale derived from 
the Zeeman signatures themselves. This is not necessarily in- 
consistent; Zeeman signatures likely trace both open and closed 
magnetic field lines, while Ca ii emission senses mostly closed 
field lines. Assuming that closed field lines concentrate around 
the equator in HD 189733 (as they do in the Sun), this would 
argue for a pole rotating slower than the equator in HD 189733. 

Given the incompleteness of our data, these conclusions are 
of course very uncertain; they will however be testable as soon 
as a data set covering the whole stellar rotation cycle is obtained. 

Folding the Ca ii H&K emission variations on half the 
planet's orbital period (ie 1 . 109 d) produces a correlation similar 
to that of Fig. |6] However, this apparent correlation is presum- 
ably the fictitious result of aliasing between the stellar rotation 
period and the length of the sidereal day. Given the sampling 
of our observations (ie roughly one series per day at fixed time, 
within about +2 hr), modulation on a period of 1 1.8 d (produc- 
ing a phase shift of -1-8.5 + 0.7% of the rotation cycle from one 
day to the next) generates a similar pattern to modulation on a 
period of 1.109 d (producing a phase shift of -10 + 7% for two 
successive observations). The phase diagram by itself is thus not 
sufficient to determine on which timescale activity is modulated. 
However, if activity was truly modulated on a 1 . 109 d period, eg 
due to tidal effects induced by the planet, we would then expect 
other activity indices (e.g. photometry traci rig surface starspots 
to show a similar behaviour; according to IWinn et aTl (12007 
and Croll et al. (2007), no modulation on a 1.109 d timescale 
is observed in their densely-sampled photometric data. We thus 
conclude that most of the activity of HD 189733 traced by the 
Ca II emission core is likely varying with the stellar rotation pe- 
riod. 

Note that we do not suggest that activity is completely inde- 
pendent of orbital phase in HD 189733; while the largest frac- 
tion of the activity level is likely driven by conventional activ- 
ity processes (and thus modulated by the star's rotation), it may 
well be that a small fraction of the observed activity is induced 
by the planet/star interaction. Given the high level of intrinsic 
variability characterising Ca ii emission (see Fig.|6]l, this option 
can only be checked from observations featuring dense temporal 
sampling over both orbital and rotational periods. 
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Fig. 5. Evolution of Ca ii H&K residual intensity in the emission 
core (in units of the unpolarised continuum), for both the 2006 
June (left panel) and August (right panel) data set. The corre- 
sponding rotation phase is indicated to the right of each profile. 



5. Transmission spectrum 

We now consider the transmission signature of the planet's 
atmosphere in the intensity spectrum, i.e., the residual on-transit 
spectrum after the off-transit spectrum has been subtracted. One 
may expect absorption signatures from the planet atmosphere 
to be present in the transmission spectrum. Similar attempts 
to detect atmospheric spectral features from the ground have 
yielded null results on the s y stem HD 209 458 (Bundv & Marc^ 
'20001; iMout ouet all 12001'; 'Br own et all 12002; Moutou et al. 
2003; Snellenl l2004r iDeming et alj |2005h . whereas space-based 
observations gave positive detections, at a level of 10"^ in the 
visible Nal doublet (Charbonneau et al., 2002) and probably 



strong er in the UV dVidal-Madiar et alj 120031 iBallester et"aLl 
120071) . 



The on-transit spectrum is the combination of all polarised 
spectra from sequences 10 and 11, i.e., in total 8 spectra of in- 
dividual SNR more than 400 at 730 nm. The final combined 
on-transit spectrum is of SNR 550 to 1100, depending on the 
wavelength. In the off-transit comparison spectrum, we included 
all spectra of the transit night except sequences 10 and 11. Each 
spectrum has to be interpolated to a common spectral format 
before combination. The impact of this interpolation should be 
negligible on features broader than the spectral sampling inter- 
val. The resulting transmission spectrum does not show absorp- 
tion features, at a level of 0.0012 relative flux in the region of the 
Nal doublet, and more generally, at a level between 0.001 and 
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Fig, 6. Equivalent width variations in the emission cores of the 
Ca n H&K lines as a function of rotation phase (derived from 
the ephemeris of eqn. |2]i, for both the 2006 June (crosses) and 
August (open squares) data sets. The equivalent width of the av- 
erage emission core is equal to 13 kms"' in HD 189733. 



0.003 over the available spectral range, from 369 nm to 1045 nm. 
Narrow Nal absorption features of 0.45% amplitude and O.Olnm 
width, such as predicted by models of HD 209458 (Tinetti et al, 
in prep) are undetected. 
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6. Planetary reflected-light spectrum 

Finally, we searched our linear polarisation spectra for signa- 
tures of light from the star reflected off the planet's atmosphere 
towards the observer. With a short semi-major axis (0.031 AU), 
a known planetary radius and the geometry of the system well 
constrained by transit observations, the system HD 189733 prob- 
ably offers one of the most favorable configurations to con- 
strain the albedo of a giant exoplanet. Upper limits on the albe- 
dos of hot Jupiters have been derived from similar searches on 
HD 75289, v And and t Boo, involvin g several tens of observing 
nights on 4m to 8m class telescopes (Charbonneau et al]|1999t 
ICollier Camero n et al. 2002; Leigh et al. 2003b a c). Their con- 
clusion is that the (assumed grey) albedos of hot Jupiters are 
apparently smaller than 0.12, i.e., much smaller than those of 
solar-system planets like Jupiter and Venus. In our case, the star- 
planet distance is the smallest among all hot Jupiters discovered 
to date orbiting stars bright enough for such observations. It thus 
maximises the strength of the expected reflected signal, making 
it larger, by at least a factor of about 2, than all systems pre- 
viously observed. The star is however fainter, compensating at 
least partly for this advantage. 

All previous attempts to detect this spectroscopic signal have 
been made using unpolarised lig ht only. However, this signal is 
expected to be linearly polarised (ISeager et al.ll200(j: IStam et akl 
[2004) and can thus also be detected through Stokes Q and U 
spectra. Detecting this signal through broad-band polarimetry 
was recently attemp ted with the new PLANETPOL polarimeter 
jHough et al.ll2006l). with which a potential signal from t Boo 
was reported ('Lu cas et"aLll2006h (at a level of about 10 ppm). 
Similar measurements can be performed through spectropo- 
larimetry. Since the reflected spectrum is expected to be po- 
larised and Doppler shifted as a result of the planet orbital mo- 



Fig. 7. Nightly averaged equivalent widths of the Ca ii H&K 
lines of the June data set as a function of rotation phase fitted by 
a cosine rotational function (top panel); residuals of the individ- 
ual Ca 11 equivalent widths after this correction is applied, as a 
function of the planet orbital phase (middle panel for the June 
data set and top panel for the August data set): the high disper- 
sion shows the activity jitter, but at such temporal sampling no 
specific correlation with the planet motion is observed. 



tion, we expect to see tiny polarisation structures in the spectrum 
of HD 189733, oscillating about the main spectral lines of the 
host star and varying in strength with orbital phase. These po- 
larisation signatures are expected to be maximum at quadrature, 
i.e., for scattering angles of about 90°. 

Although spectropolarimetric signatures of exoplanets are 
potentially easier to identify unambiguously (thanks to their 
modulated Doppler shifts) than broadband polarisation, they are 
more difficult to detect because their amplitudes are significantly 
smaller, being reduced by approximately the relative depth of 
the average line profile. Given the planet's size and orbital dis- 
tance in the particular case of HD 189733b, we expect the rel- 
ative broadband polarimetric signal to peak at about 50 ppm at 
most (assuming a grey albedo of 0.4); we thus expect the spec- 
tropolarimetric signal to reach about 15 ppm (the relative depth 
of the LSD profiles being 0.3) at orbital phases 0.75 and 0.25. 
For a grey albedo of 0.12, this signature is further reduced to 
a maximum size of only 5 ppm. Our Stokes Q and U spectra 
each show a relative noise level of 25 ppm; we thus need to ac- 
cumulate between 60 and 600 of them to reach noise levels of 
3 to 1 ppm, as is needed to detect the spectropolarimetric sig- 
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nal with a precision of at least 5cr for albedos ranging from 0.4 
down to 0.12 respectively. No useful upper limit on the albedo 
can therefore be obtained from our single pair of linear polarisa- 
tion spectra; a few tens of such observations (at the very least) 
are required to be able to derive a useful constraint. 



7. Conclusion 

HD 189733 is among the most studied planetary systems, with 
a short-period (2. 2d) Jupiter-like planet transiting the disc of 
its bright parent star. At only 0.03 1 AU distance from its rela- 
tively active host star (a K dwarf featuring moderate rotation), 
the gaseous planet is embedded in the outer corona and the mag- 
netosphere of the star. 

In order to understand the complex relationships between a 
star and a closeby Jupiter-mass planet and in particular to investi- 
gate the possible interaction that could lead to planet-induced ac- 
tivity enhancements such as those recently reported in the litera- 
ture, we decided to explore the extent and topology of the stellar 
magnetosphere of HD 189733. We used the ESPaDOnS spec- 
tropolarimeter at CFHT to look for circular polarisation Zeeman 
signatures in the line profiles of HD 189733 in 2006 June and 
August. 

Zeeman signatures are clearly detected in all spectra, demon- 
strating that a magnetic field is indeed present at the surface 
of the star. The Zeeman signatures are not modulated with the 
planet's orbital period; their temporal variation agrees better 
with the star's rotational modulation, even though the full ro- 
tation cycle of about 1 1.8 d was not fully covered by our obser- 
vations. 

Assuming that the observed variation is mainly caused by 
rotational modulation, the detected Zeeman signatures indicate 
that a large-scale magnetic field of a few tens of G is present 
at the surface of this relatively active K dwarf star. The recon- 
structed large-scale field topology is significantly more complex 
than the large scale field of the Sun, involving in particular a 
toroidal component and significant contributions from magnetic 
multipoles of order of at least 3. The magnetic topology is ap- 
parently evolving on timescales of a few rotation cycles; more- 
over, the recurrence timescale of the Zeeman signature over the 
full timescale of our observations is apparently longer than the 
estimated rotation period of 1 1.8 d. Both points suggest that sig- 
nificant differential rotation is present at the surface of the star. 

The Ca n H & K lines clearly feature core emission, whose 
intensity varies by about 10% on a day-to-day basis. This vari- 
ability is mostly modulated with the rotation cycle of 11.8 d 
and presents rapid episodes of intrinsic fluctuations. The over- 
aU modulation itself is roughly stable over the timescale of our 
observations. 

At this stage, there is very little we can state concerning 
putative activity enhancements induced by the presence of the 
giant close-in exoplanet orbiting HD 189733. In particular, its 
potential contribution to Ca ii H & K emission needs to be in- 
vestigated through further data, densely sampling both rotation 
and orbital cycles, to determine unambiguously whether, and to 
which level, the activity of HD 189733 also fluctuates with the 
planet orbital period. Ultimately, studies such as ours should al- 
low us to diagnose whether star-planet interactions rather result 
from star-planet magnetospheric friction (crucially involving the 
large-scale field of the host star) or from tidal effects (producing 
mainly comoving small-scale magnetic features with no need for 
a large-scale magnetic field on the host star); looking for corre- 
lations between the planet-induced excess activity and the pres- 



ence (and topology) of the large-scale field should for instance 
inform us on which of these 2 processes are more Ukely to occur. 
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